UNIVERSITY OF ILLINOIS 
DIGITAL COMPUTER 

LIBRARY ROUTIH E K5 - 137 
TITLE "Least Squares 

TYPE Ccanplete Program 

PURPOSE OF ROUTINE To fit a set of N weighted experimental points by the Jest 

polynomial of the form 

n-1 
F(x) = 1/2 Z a x S 
s=0 s 

The criterion of excellence for the polynomial is that the 
sum of the squares of the deviations, 



N-1 

E 

i=0 



(the f(x.) are given 



M= E [F(x.) - f(x,)ra (x.) • +-, / + x 

._ n i i w v i' experimental data) 



shall be a minimum with respect to arbitrary variations 

of the coefficients a . The resulting polynomial is known 

s 



as a "least squares polynomial". 
DURATION OF ROUTINE ,^„„2 ,__ . ,\ .. , -3 



METHOD USED 



PROCEDURE 



lONn* (n + l) x 10* ^ seconds 
The condition that M shall be minimal is equivalent to a set 
of n simultaneous linear equations for the coefficients a • 
The program produces the corresponding n x (n + l) matrix 
suitable for use with the library routines L3 or 12 (linear 
equation solver). The former should be used if the number 
of coefficients n is large, or if the equations are poorly 
conditioned . 

(1) Read in Main Tape, 

(2) Insert data prepared according to specifications 
stated belowj read it in with black switch . 

(3) To start program, throw black switch . Use of white 
switch at this point will be described below. 

(k) At end of computation and output, return to step 2 

above if it is desired to process more data; otherwise 
proceed to step 5 below. 
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(5) Clear memory, and read in the appropriate version 

of the linear equation solver. Insert output tape from first portion of program 
as data here. If more than one matrix is being solved and Routine L2 (113) is 
used, the supplement which makes L2 a closed subroutine may be used. • 

(6) The solutions to the set of n equations are, in order, 
the desired coefficients a g and a scale factor. (See description of linear 
equation solver for interpretation of the scale factor. 

SPECIFICATION OF DATA 

The data is to be prepared in three lists, the individual 
entries of which are fractions preceded by sign. 

Table I (Terminated by an "L" symbol)? The N observation 
points x . 

2*ble_TE (Terminated by an "F" symbol)? The observations 
f(x i ), listed in the same order as the x's. 

Table III (Terminated by an "N" symbol)? The weight fac- 
tors v(x ± ) in the same order as the x's. 

If all of the experimental points are assigned e<iual weights, 
Table III may be omitted. The N symbol must then be replaced by a "J". 

After the N or J is an integer specifying the number of 
coefficients n in the polynomial. This integer is terminated by spaces. Total 
space available in memory for storage is 8fk locations. Space required for 
data is 2N + 1 or 3N + 1, depending on whether or not weights are specified. 
USE OF WHITE SWITCH IN STEP 3? 

Input of data is resumed, Tables I and II being assumed 
already present in the memory. This is valuable if one has a long table of data 
which he does not wish to repreforate, but would like to alter the single 
character on the end of the tape so as to change the number of variables. At 
step 3, a short tape Jn or Nn is inserted and the white switch thrown. At stop, 
throw black switch to start program. 

If it is desired to process the same equally weighted data 
several times using different numbers of coefficients, this can be accomplished by 
reading in a short tape Jn (spaces) (n is the number of coefficients to be used 
next) following each computation, and then throwing black switch to start new 
computation. 
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SUFFLEMESTABT TBROto&IHLGBl 

If prograaraer wishes to fit the data to a function of the 

form _ -t 

o(x) « S vuW 

s«0 

(U (x) being a set of given functions), he my do so by 

s 

replacing the routine at 159 by an auxiliary. This is entered from the main 
program with s.2" 59 in A, x in 0. It is left with U g (x) in A and in 2. Scaling 
required that |u (x) | < l/2 for all x. 

S "™* 

If more than 11 words are used by this routine the storage 
locations which start at 150 must be advanced so they will not overwrite the 
auxiliary. This may be done by changing the preset parameter S6 on the master 
tape from OOF 00 150F to the desired value. 



DAT E 6/15M RT; 1/27/60 
PROGRAMMED BY R. Rubenstein 



APPROVED BY J.P.Hash 



ns 



LOCATION 


ORDER 


HOTES PAGE 1 K 3 




D.O.I 
00 3K 




3 


00 F 






00 33F 


Location of least squares routing 


k 


00 F 






00 11F 


number of decimal places in Matrix 


5 


00 F 






00 11F 


number of decimal places in solution 


6 


00 F 






00 150F 


storage 


7 


00 F 






00 QkF 


Code K2 (110 = Code Pi) 


8 


00 F 






00 139F 


Auxiliary- 


9 


00 F 






00 12F 


input routine 




26 999N 






00 33K 






Program 






00 84K 






Code N2 





00 110K 

Code PI 

00 139K 

Auxiliary 

00 12K 

Input routine 

2k 12N 



/g 
2 x -» U 



Temporary Store: 



0, 1, 2, 3, 5, 6, 7, 8, 10, 11 
-39 



Auxiliary entered with x in 0, s . 2 in A. 

Exit with x in 0, U (x) in A, 2. 

s 



I 



ummm 



3 

k 

5 

6 

7 
8 

9 

10 

11 

12 

13 

1* 

15 

16 



ORDER 

00 12K 

L5 kL 
k6 1L 
50 S6 
50 1L 
26 S7 
IA- 18L 
k2 5L 
L5 1L 
50 S6 
10 20F 
50 19L 
26 F 
k2 49S3 
50 20L 
S5 F 

26 13L 
k2 48S3 

27 12L 
L5 8S7 
kS 11L 
50 15L 
09 IF 
ho F 

L5 19S3 
IA- 8S7 
22 L 
40 1AS3 
00 kOF 
81 4F 
7^ ^1S7 
00 ^F 
91 ^F 
32 14L 
S5 F 



From 12L 



NOTES 



Plant address a 



Enter data Input 



BWSK2 



13 



Switch after data is read in 

Leaves first location of previously 
read in data in A. 



w end of weights 

a 

put order in Q. 
store Q at 1*«33 
end data 



end of f , 



L end of x. 

Store -1 after last x J 



Plant address of next data storage, 



read in constant 



LOCATION 
IT 

18 

19 
20 



ORDER 

42 50S3 
24 S3 

26 1L 
00 6L 

27 18S3 
75 F 
40 2F 
75 F 



to main 

program 



HOTBS 



store n. 



nm 3 



K5 



use without weight table 



use with weight table 



H-l 



Punches out matrix 



A rs "£ V* V x i>V x l> 



r = 0,l 1 2,...,n-l 
s = 0,l,2,...,n-l 
r = 0,1,..., n-1 

Each row is scaled to < l/2 but greater than ~ l/8 



3 
4 

5 
6 

7 

8 



92 993F 
27 2L 
92 770F 
F5 6F 
42 6F 
LO 50L 
36 6L 
*U 8F 
L5 26L 
42 18L 
42 19L 

23 33L 
L7 48L 
00 16F 
82 4F 

92 1021F 

24 S9 



from 
29L 
From OL 



From 3L 



Print spaces 

r -»0 

Print N (End of row) 

r -> r + 1 



r-n 

< Clear store. 

Plant address of shift orders 

s ->0 

Print character at end of matrix 



Spaces 

On to next data 



t LOCATION 


ORDER 




NOTES PAGE 4- 




LI 50L 


From 31L 




9 


Ik JF 








32 4^L 




s - n 


10 


L5 TF 
50 10L 




< 


11 


2o So 




Calculate U • 

B 




40 3F 




Store U at 3. 
g 


12 


15 .6f 

50 .121, 






13 


26 SS 
50 3F 




Calculate tJ 
r 


14 


40 2F 




Store U at 2 
r 




75 F 


By JL, 


• v 


15 


40 10F 
S5 S6 


23L 


Calculate U • w 
s 


16 


■40 11F 
50 UF 




Double Precision 


IT 


7J 2? 
50 10F 






18 


Jk 2F 








10 F 


By 4L, 


ITU, . 2~ P . 2 U . w 
r s 


1 9 ,. 


- 40 IF 


5L # 25L, 






00 F 


43L 




20 


14 5F 




Store accumulated sum 


21 


ho 5F 
F5 27L 




Eu^jju^).^ -2^ .2 tt ' .. 


- 


42 27L : 




Prepare to treat next experimental 


> ■$$'^ 


- : F5 Hi. 
.42 44l 




points* 


23 


F5 14L 
42 14-L 






24 


IL 5F 
32 27L 




Test scaling of matrix element 


25 


F5 18L 




Scale down on next run 




42 181 








K3 



I ■ : — ; — 

LOCATION 


ORDER 


NOTES PACE 5 


2.6 


L5 5F 


By 21L, 35q 




10 IF 






27 


ho 5F 


From 38L 






L5 F 




Bring down next x. -* 


s8 


kO F 


From 24L 


*■ m 




L7 8F 




Test for maximum of terms in above 


29 


L2 5F 
36 31L 




sum. 


30 


L5 5F 
40 QF 




Store new maximum value 


31 


L7 F 
32 8L 


From 29L 


Test for end of data 


32 


LI 4-8L 




Switbn to print every other 




32 4-5L 




time 


33 


F5 7F 








4-2 7F 


From 5L 


s -» s +1 


3h 


FO 50L 








32 38L 




s - n 


35 


L5 15L 




< 




4-2 2 7L 




Store address first x 

1 


36 


L7 4-8L 




Store address first f 




k-2 kkL 




37 


L5 ^9L 




Store address first w 




42 14L 




38 


27 27L 




N(5) ->0 ' 




Li 4-8L 


From 3^L 


Try switch again, and 


39 


40 48L 








36 1L 




change ito 


4o 


L3 8F I 




Scale 




12 &F 


From kkL 


up every other time through 


in 


32 5L 
00 IF 




the row of matrix 


4? 


32 5L 
40 8F 






4-3 


F5 19L 


J Scale up on next run "by factor 2 • 


k- L 


_ 


• 



K 3 



LOCATION 

hi 

k8 

50 



ORDER 
'te'19L 
23 ^OL 
LI F 
22 11L 

L5 5? 
50 S^ 
50 i^6L 
26 26S7 
26 33L 
00 S5 
00 F 
00 S6 
00 F 
00 F 
00 F 
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NOTES 



By 22L, 36l| N(A) =0 
From 9L -f goes to A 



From 32L Print 

A 

rs 

back for next matrix element 



storage of constants 



(number of variables) 



1 



K 3 



1 



